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METHODS AND SYSTEMS FOR IMPROVED EXTRACTION OF A SIGNAL 

FROM A NOISY ENVIRONMENT 

Field of the Invention 

The present invention relates generally to signal processing methods and systems, more 
particularly, to methods and systems for improving the signal-to-noise plus interference ratio of 
a signal in a noise/interference backgroimd. 
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Background 

The problem of extracting a usable signal from a noisy environment is a common one in 
the field of signal processing. Noise is generally defined as an unwanted or interfering compo- 
nent that degrades system performance and detracts from the processing of the desired signal. 
Noise may be provided by numerous sources both external and intemal to a signal processing 
device or receiver. In some applications, an incoming signal may be corrupted by strong additive 
noise such as, for example, weak radar echoes from a target in the presence of strong back- 
ground clutter. In some cases, signals traveling over existing commimications systems experi- 
ence phase and amplitude distortions from channel noise that is intemal to the receiving device, 
much like tinnitus in the human ear or thermal noise in an optical or radio frequency (RF) detec- 
tor. In still other applications, noise may be introduced into a received signal by deliberate and 
intentional interference, such as in jamming. Jamming is any intentional or deliberate signal 
interference utilized by hostile forces to disrupt normal receiver operation. 

Because of the noise presence, receiving systems generally exhibit a sensitivity threshold 
that limits the processing of signals to those which exceed the noise level. The utilization of 
information contained in signals that are obscured or contaminated by noise is limited by the 
ability of the receiver processing device to separate the signal components from the noise, elim- 
inate the noise component, or enhance the signal to a level in excess of the noise to permit the 
extraction of the desired information. 

Radar performance, in particular, may be severely degraded in the presence of noise, 
"clutter," jamming or chaff Clutter is a term used to describe the return emanating from objects 
which tend to obscure the signal associated with the "true" target return. Clutter can be classified 
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in two main categories: surface clutter and airborne or volume clutter. Surface clutter includes 
the radar retum from trees, vegetation, ground terrain, man-made structures and the sea surface. 
Volume clutter may be produced by large ensembles of airborne scatterers such as that associ- 
ated with chaff, rain, birds, and insects. Chaff consists of a large number of small dipole reflec- 
tors that have large radar cross section (RCS) values. It may be released by hostile aircrafts or 
missiles as a means of an electronic countermeasure (ECM) in an attempt to confuse a radar. 
Clutter is particularly troublesome for most radar detection systems because it introduces a 
noise-like retum with a random phase and amplitude that is may exceed that of the received sig- 
nal. Unlike clutter, so-called "white noise" generally introduces the same amount of noise power 
density across the entire operating bandwidth of the radar. 

There are several known techniques for extracting signal information from noisy envi- 
ronments caused by one or more of the situations described above. Known techniques include, 
for example, matched filtering, KahnanAVeiner filters, moving target indicator (MTI) systems, 
and coherent integration. Each of these techniques, while suitable for some applications, suffers 
from limitations, as described in more detail below. All of these techniques are described in 
more detail in Bassem R. Mahafza, "hitroduction to Radar Analysis" (CRC Press, 1998) and in 
Bassem R. Mahafza, "Radar Systems Analysis and Design Using MATLAB" (CRC Press, 
2000). 

In a system which utilizes a matched filter, the signal-to-noise ratio (SNR) of the 
receiver output is enhanced by or through correlation of the received signal with the transmitted 
wave form. In a radar system, for example, a received signal consists of an amplitude weighted 
and time delayed version of the transmitted signal plus noise. The receiver processes a replica of 
the transmitted signal together with the received signal to determine the correlation between the 
two signals. Maximum correlation indicates the location of the desired signal, i.e. the "true" tar- 
get retum, within the receiver time window. While matched filtering may work in many situa- 
tions, it may be inadequate in clutter, chaff, and jamming environments. In these situations, the 
amount of interference signal is so overwhelming that processing requires prohibitive pulse 
lengths or system band widths to recover the true signal from the background interference. 

Kalman filters use Markov processes to estimate the desired signals and, as such, are dif- 
ficult to implement because of the mathematics involved. Furthermore, the proper selection of 
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gain parameters requires accurate a priori knowledge of the target dynamics and is not straight- 
forward. Weiner filters, a subclass of Kalman filters, are effective with thermal background 
noise but are not generally effective against clutter or jammer retxims. 

In coherent integration systems, otherwise known as pulse integration systems, coherent 
radar signals are added together to increase the SNR. Two or more radar signals are said to be 
coherent if the amplitudes and relative phases of the signals have a known relationship even 
though they may be separated in time. If the noise is uncorrelated, the SNR increases linearly 
with the number of signals integrated by the processor. The accuracy of these systems are lim- 
ited by their ability to integrate a sufficient number of samples over an extended period of time 
since fluctuations in the amplitude and phase of the target return may cause the desired signal to 
decorrelate and become noncoherent. 

In moving target indicator (MTI) systems, delay-line cancelers are used to suppress tar- 
get-like return signals produced by clutter, and allow return signals from desired targets to pass 
through with little degradation. Clutter is distinguishable fi^om receiver noise by its relatively 
low-fi:equency, narrow spectrum. An MTI system cancels the clutter retimi fi"om the incoming 
signal by subtracting successive echoes fi-om the same location. The problem with MTI systems 
is that they exhibit "blind speeds," corresponding to Doppler firequencies for which the radar is 
unable to detect targets. 

Fig. 1 depicts the output gain of a three-pulse MTI canceler along the vertical axis with 
frequency shown along the horizontal axis. As shown in Fig. 1, blind speeds occur when the tar- 
get velocity results in a Doppler fi-equency equal to an integer multiple of the radar pulse repeti- 
tion frequency which, in this example, is at 100 Hz, 200 Hz, 300 Hz, etc. A target, having a 
relative velocity corresponding to a Doppler frequency equal to an integer multiple of the pulse 
repetition frequency, will not be detected. Another Umitation of MTI systems is that the pass 
band response (gain) may vary widely over the bandwidth of interest which masks the effects of 
a change in the target velocity. 

The performance of MTI systems can be improved by increasing the number of delay 
lines. Fig. 2 shows an example of a MTI system with multiple delay- line cancelers and a 
weighted summer. As shown in Fig. 2, the signal is tapped at multiple locations at even intervals 
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delay time, T, seconds apart. The delayed pulses are weighted by binomial or other coefficients 
and then added together to produce the output of the MTI filter. 

The perfomiance of this type of MTI system is known to be dependent on the choice of 
weighting coefficients. There are numerous methods for choosing the weights, some of which 
are discussed in Fred E. Nathanson, "Radar Design Principles," (McGraw-Hill Inc., 1991, 2nd 
ed.), pp. 410-12. One common weighting scheme is to choose weights that are binomial coeffi- 
cients (coefficients of the expansion (1 + a)'') with altemating signs. Notwithstanding the 
effectiveness of MTI filters with binomial coefficients, these systems still exhibit blind speeds 
and wide pass-band (gain) fluctuations may mask the effects of changes in the target velocity. 
Therefore, research continues to define alternate weighting schemes and associated coefficients 
to minimize the effect of blind speeds and provide a better pass-band response. 

Summary of the Present Invention 

Methods and systems consistent with the present invention allow identification of a true 
signal contained in a signal containing the true signal and noise. In general, digital signal infor- 
mation representing a signal of interest plus noise is utilized by the present invention. The first 
A'^ samples of digital signal information are stored with the Nth sample being stored in manner 
which renders it accessible for additional operations. A specially selected set of weights are 
applied to the buffered N samples and, additionally, phase rotation is applied to the Nth sam- 
ple. The phase rotated Nth sample and weighted samples are combined using a first equation, 
described in more detail below. The resulting signal, which exhibits an increased Signal-to- 
Noise ratio (SNR) and may be more effectively utilized in subsequent MTI processing by virtue 
of the operations performed on the previous N samples as described herein, is then available for 
further processing using conventional techniques. 

Brief Description of the Drawings 

The accompanying drawings, which are incorporated in and constitute a part of this 
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Fig. 1 is a graph showing the output of a three-pulse MTI canceler indicating blind 

speeds; 

Fig. 2 is a block diagram of a MTI system with multiple delay-line cancelers and a 
weighted summer; 

Fig. 3 is a flow chart illustrating a method for improving signal-to-noise ratio (SNR) in a 
signal containing noise consistent with the present invention; 

Fig. 4 shows one exemplary method for obtaining samples consistent with the present 
invention; 

Fig. 5 is a comparative graph showing blind speeds experienced by both a prior art sys- 
tem and the present invention; and 

Fig. 6 is a schematic block diagram of a signal processing system consistent with the 
present invention. 
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Detailed Description 
A. Overview 

Methods and systems consistent with the present invention allow detection of a signal 
contained in a noise environment. In general, a received analog signal containing information of 
interest immersed in a noise/interference background is coherently processed to resolve it into 
its complex quadrature components. The quadrature components of the analog signal are sam- 
pled at predetermined time intervals, , and, through an analog-to-digital converter, expressed 
as a digital stream of the sampled signal quadrature components. consecutive samples of the 
digitized signal quadratxire components are stored. A specially selected set of weighting coeffi- 
cients are applied to the samples. Additionally, phase rotation is applied to the last, i.e.. Nth 
sample. The phase-rotated Nth sample and the N weighted samples are combined using Equa- 
tions 1 through 4 described in more detail below. The resulting modified signal may then be 
processed by conventional signal processing means. A next set of A'^ samples are then processed 
as described above to yield modified signals every A/ . 
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The present invention is an improvement over prior art systems previously mentioned 
because, while simple to implement, it can be easily integrated into existing systems. The 
present invention allows for fast and accurate detection of signals that other systems can not 
detect or which may only be able to detect after a lengthier processing period or by using a com- 
plex processing scheme. Furthermore, methods consistent with the present invention may be 
used to improve the detection of signals in a broad range of applications, including but not lim- 
ited to radar systems and other communication systems. Methods consistent with the present 
invention may be used to improve signal detection in optical and sonar devices. 

Reference will now be made in detail to exemplary embodiments of the invention that 
are also illustrated in the accompanying drawings. The present invention may be implemented 
using hardware, software or a combination of hardware and software. The description includes 
exemplary embodiments, other embodiments are possible, and changes may be made to the 
embodiments described without departing from the spirit and scope of the invention. Wherever 
possible, the same reference numbers will be used throughout the drawings to refer to the same 
or like parts. 

B. Method 

One embodiment of the present invention consists of a radar processing method which 
provides significant enhancement of the SNR at the receiver device through a unique type of 
coherent processing. More specifically, one embodiment of the present invention comprises per- 
forming a special weighted averaging of the quadrature components for a number of pulses. The 
method of the present invention produces new estimates of the pulses removing from the signals 
a significant portion of the spatially incoherent components emanating from noise, clutter, or 
jamming sources. 

Fig. 3 is a flow diagram of a method for detecting a desired signal in a signal containing 
the desired signal and noise consistent with the present invention. In a radar appUcation, the ana- 
log radar signal is received (step 305). The received radar signal may be expressed in the form 

x(t) = A(t)e^^^^ ^^^^^ , where ©q is the intermediate carrier frequency (IF) of the radar and 
the A(t) and (p(t) represent the amplitude and phase modulation of the target return, respec- 



6 



Atty, Dkt. No. : 7607.0001-00 



tively. The received radar signal is coherently processed to resolve it into its in-phase Xj{t) and 
quadrature components given by x(t) = Xj(t) + ixgit) . After coherent processing, the 
signal is converted to a digital format (step 310). 

In step 3 1 5, the window size, (A^- 1 ) A/ is determined. Window size is the time period 
over which the samples will be taken. Window size is a function of the application and is based 
on, for example, a sampling rate and the characteristics of the signal. For example, a smaller 
number of samples or a faster sampling rate may be required for a faster moving target to ensure 
that the target retum does not completely decorrelate within the window. Initially, a window 
size may be set based on anticipated target dynamics or the known characteristics of a target of 
interest, however, the window size may also be adjusted manually or dynamically based on the 
characteristics of the incoming signal. 

Fig. 4 shows one exemplary method for obtaining signal samples consistent with the 
present invention. The upper graph plots the in-phase component of the analog signal. The 
lower graph plots the quadrature component. In the example shown in Fig. 4, four samples are 
taken (A^ = 4) every At time period. The total window size therefore is (A'^- 1)A^ . The sam- 
pling interval corresponds to the system's A/D sampling rate, which depends on the charac- 
teristics of the signal of interest. As stated above, the samples shown to the left of the graphs, 
Xj^ - to- Xj^ and Xg^-to- Xq^ are the in-phase and quadratxxre components of the analog sig- 
nal, ;c(/) . Each of these samples may then be quantized, thereby converting the analog signal to 
digital format. 

The N digital samples of the analog signal may be stored, for example, in memory (in a 
software implementation) or a shift register (in a hardware implementation) (step 320). Each 
sample consists of a pair of digital quadrature component values Xj and Xq , as shown in Fig. 4. 

IfN- 4 , then complex samples Xj, x^ and x^ are stored. Those skilled in the art will rec- 
ognize that the in-phase and quadrature components may be stored separately or together as x(t). 
Additionally, particularly in a hardware implementation, the last sample, x^, which in this 

example is x^ , may be stored in a second buffer for further processing, (step 325). 
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The last sample, Xj^, is processed using phase rotation (step 330). One type of phase 

rotation is described in Colher, Robert J., et al., "Optical Holography", (Academic Press, 1971), 
pp. 22-24. Phase rotation may be performed by, for example, calculating the conjugate in soft- 
ware or performing an equivalent bit manipulation in hardware. 

In step 335, specially selected weighting coefficients, , are applied to the stored sam- 
ples. Each of the weighting coefficients belongs to a special class of integers that possess imique 
characteristics. Use of the unique relationships between these weights allows for maximizing the 
ratio of the desired signal power to the undesired signal power, that is, maximizing the SNR. In 
one embodiment of the present invention, the weighting coefficients are comprised of members 
of the class of integers described as Mersenne primes, or their reciprocals. A Mersenne prime is 

a prime niunber, , where = 2" - 1 and n is also a prime number. The set of Mersenne 
primes, {MJ , includes {3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 

3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 
216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377,...}. In one embodiment of the 
present invention, four consecutive Mersenne prime numbers are used. 

The weights are applied by modifying the buffered signal samples with the weighting 
coefficients by mixing or multiplying. 

In step 340, the processed last sample, jc^, is combined with all weighted samples to pro- 
duce a new sample estimate, that is, modified signal, j;^,. In methods consistent with the present 

invention, the phase rotated Nth sample and the N weighted samples are combined. A new esti- 
mates of the incoming signal samples is formed and declared by 3;, as described below for 
N = 4, The first three sample estimates are: 
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= + "3 = ^3 



Eq(l) 
Eq(2) 
Eq(3) 
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where s is the desired signal return and u is additive uncorrelated noise with zero mean and 
variance (noise power). Note that the desired signal s and the noise o are also spatially inco- 
herent. More precisely, ^[^,Oy*] = 0 = £'[u^Uy*] , where ] refers to the expected value 

operator and the asterisk denotes complex conjugation. Thereafter, for / = A^> 4 , the new 
samples y^; i>N = 4 arc computed as follows 



X; 



Ml M2 M3 M4 



Eq(4) 



where M, > M2 > M3 > are weights as discussed further above. The term x^* is the com- 
plex conjugate of x . . In some embodiments of the present invention, the coefficients in Equation 
4 take on alternating signs. 

Denote the new complex sample estimates y. as 



Eq(5) 



It follows that the following is true 



which is the desired signal power. And that 



Eq(6) 



E[s,vj*] = E[hvj*] = £[5,u/] = Eihv*] = 0 



^[u,u/] = £[u,o/] = 0 



for all (ij) 



Eq(7) 
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£[u,u/] = ^[u,o/] = E[v,Vj ] = E[ViVj ] = 0 for all (i^j) 

^ for (i=j) 



Eq(8) 



After combination, the resulting modified signal is output (step 345). The outputted signal may 
be stored or transmitted to a data processor for further processing of the signal using conven- 
tional techniques. 
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If the process detects the end of the signal stream (step 350), the process terminates. If 
there are more signal samples remaining to process, the buffers holding the signal samples are 
shifted (step 355). The last sample, , is shifted out of the shift register and a new sample is 
stored. As long as there are more samples to be processed, in step 355, the buffer is shifted, and 
the process proceeds with processing the next set of samples. 

The present invention significantly enhances the SNR over conventional methods. The 
improvement in the SNR can be determined from E[y^] . More precisely (for N = 4), 

where (w^ = l/Mf), for I = 1, 2, 3, 4 . Substituting = + and Equation 5 into Equa- 
tion 9 yield 



E[yi] = E[(s,^V,r{ (mj(5,- + 0,-) + '«2(^/-l+0,-l) + ^3(^/-2 + ^i-2) + 

/«4(^r-3 + ^/-3)}] 

Utilizing Equations 6, 7, and 8 into Equation 10 gives 



Eq(lO) 



2 2 2 

E[y-] = m^\s\ +m^G^ +H (^2 + ^3 + ^4) 

It follows that the SNR is now computed as 



Eq(ll) 



(SNR),s = 10 X log[-J L_J 



Eq(12) 
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The ratio defined in Equation 12 is maximized when (m j + ^2 + /M3 + m^) »m^. This condition 
is attained by selecting the weights M3, as prime numbers. This ratio is fiirther 

maximized when these weights are chosen firom the class of prime numbers described in the 
foregoing. 

An additional feature of the present invention relates to the improvement of MTI pro- 
cessing. Fig. 5 is a comparative graph showing blind speeds experienced by an MTI system and 
the present invention. As shown in Fig. 5, the present invention results in significant improve- 
ment in normalized gain, which significantly improves signal detection capabilities. 
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C. System 

Fig. 6 is a schematic block diagram of a signal processing system consistent with the 
present invention. As shown in Fig. 6, one embodiment of the present system includes an 
antenna 610, a receiver 620, an A/D converter 630, a signal processor 640, and a data processor 
650. 
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Antenna 610 collects a modulated signal, such as a radio frequency signal, and passes it 
to receiver 620. The initial stages of Receiver 620 may comprise a low noise ampUfier, an RF 
mixer, an intermediate frequency (IF) filter, bandpass filter, mixer, low pass filter, optional gain 
filters, and other equipment capable of changing the signal to a baseband signal. The signal is 
then coherently processed to resolve it into I and Q components. The signal is then passed to PJ 
D converter 630. 

In a system consistent with the present invention, A/D converter 630 is comprised of 
analog-to-digital processing circuitry and filtering components, as is well understood in the art. 
A/D converter 630 converts the analog signal into a digital signal at a suitable sampling rate. 
Any sampling rates may be used provided that appropriate digital signal processing for the par- 
ticular application is performed. For example, a fast moving radar target may require a higher 
data rate. For speech, you need only twice the bandwidth of human speech. A/D converter 630 
may, in fact, be any standard analog-to-digital converter means capable of receiving an analog 
electrical signal and outputting a digital signal consistent with the selected sample rate. 

Signal processor 640 may be a standard or specialized computer, depending on the appli- 
cation, implementing the method of the present invention. Signal processor 640 may be any 
commonly available digital processor such as, for example, a Pentium. Signal processor 640 
may be a single processor or multiple processors. As described above, the method of the present 
invention may be implemented in software, hardware, or a combination thereof 

Data processor 650 is a standard computer or special hardware. Data processor 650 
could be any commercial or standard computing device, or a specially designed computing 
device that does decision processing. Data processor 650 may be any commonly available digi- 
tal processor such as, for example, a Pentium. Data processor 650 may be a single processor or 
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multiple processors. In a radar system, for example, the system would process the resulting sig- 
nal to detect the presence of the target, and determine its range, velocity, and other target param- 
eters. 

The foregoing description of an implementation of the invention has been presented for 
purposes of illustration and description. It is not exhaustive and does not limit the invention to 
the precise form disclosed. Modifications and variations are possible in light of the above teach- 
ings or may be acquired from practicing of the invention. For example, the described implemen- 
tation includes software but the present invention may be implemented as a combination of 
hardware and software or hardware alone. The scope of the invention is defined by the claims 
and their equivalents. 
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